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ABSTRACT 

Aims. Estimating Alfven speeds is of interest in modelling the solar corona, studying the coronal heating problem and understanding the 
initiation and propagation of coronal mass ejections (CMEs). 

Methods. We assume here that the corona is in a magnetohydrostatic equilibrium and that, because of the low plasma ji, one may decouple the 
magnetic forces from pressure and gravity. The magnetic field is then described by a force-free field for which we perform a statistical study 
of the magnetic field strength with height for four different active regions. The plasma along each field line is assumed to be in a hydrostatic 
equilibrium. As a first approximation, the coronal plasma is assumed to be isothermal with a constant or varying gravity with height. We study a 
bipolar magnetic field with a ring distribution of currents, and apply this method to four active regions associated with different eruptive events. 
Results. By studying the global properties of the magnetic field strength above active regions, we conclude that (i) most of the magnetic 
flux is localized within 50 Mm of the photosphere, (ii) most of the energy is stored in the corona below 150 Mm, (iii) most of the magnetic 
field strength decays with height for a nonlinear force-free field slower than for a potential field. The Alfven speed values in an isothermal 
atmosphere can vary by two orders of magnitude (up to 100000 km/s). The global properties of the Alfven speed are sensitive to the nature of 
the magnetic configuration. For an active region with highly twisted flux tubes, the Alfven speed is significantly increased at the typical height 
of the twisted flux bundles; in flaring regions, the average Alfven speeds are above 5000 km-s~' and depart strongly from potential field values. 
Conclusions. We discuss the implications of this model for the reconnection rate and inflow speed, the coronal plasma [} and the Alfven transit 
time. 
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1. Introduction 



In|1942, Alfven discovered that magnetohydrodynamic (MHD) 
waves are of three kinds: two magnetoacoustic waves (slow 
and fast modes) which are compressible and can be subject 
to damping, and the so-called Alfven wave which propagates 
along the magnetic field and is incompressible. The Alfven 
wave propagates at the Alfven speed given by: 



temperature of the corona (about 1 MK) compared with the 
photosphere (about 4600 K). It is thought that the heating 
mechanism for the corona is of magnetic origin: the magnetic 
energy stored in the corona is released and converted into heat 
either by magnetic reconne ction or by dissipation of waves (see 
review bylfQimchuk 2006, and reference therein). The Alfven 
speed appears as a tuning parameter in both heating processes. 
For instance, in a reconnection model such as the Sweet-Parker 



B 



VA 



(1) 



where B is the magnetic field strength and p is the density of the 
plasma. Since Alfven's work on MHD waves, the Alfven speed 
has been shown to be an important quantity in understanding 
many physical processes in the solar atmosphere such as (i) 
the coronal heating problem, (ii) the (non-) equilibrium state of 
the corona, (iii) the initiation and propagation of coronal mass 
ejections (CMEs). 

The coronal heating problem aims basically to identify and 
understand the mechanisms responsible for sustaining the high 
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model (Sweet 1958; Parker 1963), the outflow speed is of the 
order the Alfven speed of the inflow region. Therefore, the 
measurement of the local Alfven speed in the corona is an es- 
timate of the outflow speed in t he case of a possible recon- 
nec tion process. Observationallv. | Narukage & Shibata (2006) 
and iNagashima & Yokovamal ( 120061) have analysed numerous 
EUV and soft X-ray flares to deduce an inflow speed of a 
few tens of km-s and a reconnection rate between 10~ 3 and 
10 _1 of the Alfven s peed. These results are consistent with 
the estimates made by iDere ( 19961) . In terms of MHD waves, 
the mode conversion between slow and fast magnetoacoustic 
waves is the efficient when the sound speed is equal to the 
Alfven speed for a disturbance propagating along the mag- 
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netic field. Nevertheless under coronal conditions, the heating 
by dissipation of Alfven wave is more likely than by mag- 
netoacoustic waves. The dissipation of shear Alfven waves 
generated in the low corona can be considered as a plausible 
heating mechanism especially when the d ynamics of Alfven 
waves is dominated by phase-mixing (jH e vvaerts & Priest 
ll983llHood et al.ll997[lNakariakov et al.ll997h or resonant al> 
sorption dGrossmann & Smith! 1 9881: IPoedts et al.ll 1 989L 1 1 990h . 
The damping of shear Alfven waves by phase-mixing is asso- 
ciated with the existence of a gradient in the Alfven speed and 
inhomogeneities in the background stratification. Resonant ab- 
sorption is an other viable mechanism for coronal heating when 
the driving velocity of the magnetic structure is equal to the lo- 
cal Alfven speed. Recently iMcLau ghlin &Hoodl (12004 l2006h 
showed that fast magnetoacoustic waves can be dissipated near 
a null point in a low-/3 plasma and that Alfven waves are more 
likely dissipated near separatrices. Both fast and Alfven waves 
contribute to Joule heating in the corona localized nea r topolog- 
ical elements. Furthermore. lLongcope & Priest! (120071) showed 
how fast magnetoacoustic waves maybe launched by transient 
reconnection. 

A coronal magnetic configuration can be considered to be 
in equilibrium if the time scale of the evolution r eq is longer 
than both the Alfven time scale ta (or Alfven transit time) and 
the reconnection time scale T rec (e.g., Antiochoslll987t IPriestl 



1999). For a typical length of 100 Mm and an Alfven speed 
of 2000 km-s , ta is of 1 min. The reconnection time r rec is 
of few seconds. The transit time T eq is often assumed to be 
of the order of 10-20 min in the corona for a typical length 
of 100 Mm and an upper limit of the plasma velocity of 100 
km-s -1 corresponding to a sub-sonic coronal flow. In the limit 
of T eq » ta » r rec , the magnetic configuration can relax 
to a mini mum energy s tate, i.e. a linear force-free field ac- 
cording to lWoltjerl (119581) . Then the evolution of the magnetic 
field can be described as a series of linear force-free equilibria 
(Hev vaerts & Priest! 1 1 9841) . If T eq > ta » T rec , the magnetic 
configuration can relax to a nonlinear force-free equilibrium. 
Reg nier & Can field (2006) showed that the evolution of an ac- 



tive region can be well described by a series of nonlinear force- 
free equilibria, the time span between two equilibria being of 
15 min. 

In the different classes of CMEs, two main ingredients are 
often considered: the existence of a twisted flux tube or sheared 
arcade that can store mass and magnetic energy, and the exis- 
tence of a current sheet above or below the twisted flux tube 
or sheared arcade in order to change the connectivity of the 
field lines by magnetic reco nnection. Different models exist . 
The classical CSHKP model JCarmichaelll 964 Isturrockl 1968 ; 
Hirayamal 1974 ; Kopp & Pneuman 1976h assumes the forma- 



tion of a current sheet below a rising twisted flux tube in the 
pre-flare phase leading to the eruption of th e flux rope and 
the fo rmation of post-fla re loops. According to Cargill & Priest 
(1982) and contrary to iKopp & Pneuma nf (1 19761) . the evolu- 
tion of post-flare loops is better describe d by slow-shock waves 
(ohmic heating). In th e breakout model ( Antiochos et alJI 1999: 



Mac Neice et alj|2004l) . the increase of shear in the underlying 



field lines of a quadrupolar configuration leads to field lines 
opening at the top of the magnetic system. This model leads 



to the formation of a flux rope in 3D ( ILvnch et al.ll2005l) . The 
development of sheared arcades and/or the instability of a flux 
rope are then a prerequisite to a CME. The coronal sheared ar- 
cades are formed by shearing motions of magnetic polarities in 
the photosphere. The formation of unstable flux ropes and/or 
the destabilisation of an existing flux r ope in the corona ar e 
often considered: kink instability (e.g . . iHo od & Priest 1981), 



flux emergence (e.g.. lAmari et al.l2000tlChen & Shibata 2000), 



cancellation of flux (e.g., Parnell et al. 1994; Priest et al 



1994), 



tether-cutting (e.g jMoore & Labo nte 1980). As noticed above, 
the reconnection process is an important step in the CME: the 
speed and the acceleration of a CME is related to the local 
Alfven speed and to the reconnection rate. Therefore, the na- 
ture of CMEs can be inferred from the determination of the 
local Alfven speed: to know where the reconnection is more 
likely to happen and to estimate the speed of the CME. 

To our knowledge, few attempts have been made to char- 
acterise the Alfven speed in a coronal magnetic configuration. 



Derd (119961) performed a study of the Alfven speed in active 



regions, quiet-Sun regions and coronal holes based on mag- 
netic field estimates and/or potential extrapolations, the den- 
sity and loop length bein g derived from coronal observations. 



Warmuth & Mann d2005l) derived the Alfven speed associated 



with the propagation of wave disturbances by modelling the 
magnetic field of an active region by a magnetic dipole su- 
perimposed on that of the quiet Sun as in iMann et aL I d2003l) 
and c onstraining the density by observations. Following iLin 
(2002), the local Alfven speed in an isothermal atmosphere 
with a constant gravity is assumed to decrease with height un- 
til about 1.5 solar radii (~ 350 Mm above the surface) and 
then to increase. When a solar wind component is add e d to 
the m odelled atmosphere (ISittler & Guhathakurtalll999l) . ILin 
d2002l) showed that the Alfven speed decreases with height but 
is still consistent with Alfven speed values in the low corona 
(< 350 Mm). Nevertheless, the variation with height and the 
dynamical range of the local Alfven speed are not well-known 
in the low corona. 

As the physical processes mentioned above rely on the na- 
ture of magnetic configurations in the corona, the first part of 
this paper is dedicated to a determination of the magnetic field 
strength and its variation from the photosphere to the corona. 
The coronal magnetic field is not easily measurable. Several at- 
tempts have been made to me asure the magnetic field in promi- 



nences via the Hanle effect dLeroy et al.lll98 4). or above ac 



tive r egions from gyro-emission frequencies dWhite & Kundu 



19971) . Recently the possibility of measuring the magnetic field 



based on infrared wavelengths has been investigated as well 
as the inference of the corona l magnetic field s trength from 



the pro perties of type II bursts (ICho et al.ll2007l) . iPoletto et al 



(119751) derived the magnetic field strength of several active re- 
gion field lines using a potential field extrapolation. Based on 
this method the magnetic field strength decreases with height 
faster than 1/r from the photosphere and depends on the to- 
tal magnetic flux and the spatial distribution of polarities on 
the photosphere. The significant results in measuring the coro - 
nal magnetic field have been summarized by Koutchmy ( 2004 ). 
Using a nonlinear force-free field method to extrapolate the ob- 
served photospheric magnetic field into the corona, we study 
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the global properties of the magnetic field in the corona for dif- 
ferent active regions. In lRegnier & Priest d2007allbh . the authors 
described the geometry of field lines as well as the variations of 
energy density with height. In this paper, we focus on the dis- 
tributions of the magnetic field strength with height revealing 
different behaviours between different active regions. 

The second part of this paper is to derive Alfven speeds 
in the corona. The assumption of our model is that the corona 
satisfies a magnetohydrostatic equilibrium given by: 



Vp+pg + j AB = 0, 



(2) 



describing the balance between plasma pressure p, gravity 
and magnetic forces. Several magnetohydrostatic models al- 
ready exist and are based on the extrapolation of the mag- 
netic fie ld into the c orona satisfying Eqn.Q with or without 
gravity Cowl 119851 iBogdan & Lowlll98fAlLow|[l99ll Il992 . 
1 99 31 lNeukirchlll995l ILowI 120051: lNeukirchlll997l). Recently 
Wiegelm ann & Neukirch ( 20061) have developed an optimiza- 
tion scheme to address this problem. If we assume that the mag- 
netic forces dominate the corona, so that there is a force-free 
equilibrium for the coronal magnetic field and a hydrostatic 
equilibrium for the distribution of the plasma along each field 
line: 



(V A B) A B = o(fi) 
and 

B-(-Vp+pg) = 0. 



(3) 



(4) 



These equations imply that there is no feedback of the coro- 
nal plasma on the magnetic field. However the magnetic field 
strongly influences the plasma by determining the shapes of 
coronal loops. Under coronal conditions, we assume that the 
plasma /3 is much less than 1 : Eqn (01 is a force-free equilibrium 
and Eqn. (0]i is the hydrostatic equilibrium. Solving these equa- 
tions, we investigate the values and distribution of the Alfven 
speed in the low corona assuming that the corona is isothermal. 
We only consider the global properties of the Alfven speed and 
not the Alfven speed profiles along individual field lines. 

In Section |2] we first describe the observed active regions, 
the magnetic field model as well as the isothermal atmosphere 
model required to compute the Alfven speeds. In Section [3] 
we then study the field strength and Alfven speed in a twisted 
bipolar field. We apply this model to four solar active regions 
in Section|4] We discuss the implications of these results on the 
physics of active regions and CMEs in Section[5] 

2. Observations and models 

2.1. Observed active regions 

To show the diverse and complex nature of the coronal mag- 
netic field, we have selected four different active regions with 
different types of activity (confined flares, flares associated 
with a CME or filament eruptions) and at a different stage of 
their evolution (before or after a flare): 



flux and magnetic polarities diffusing away). A filament 
eruption associated with an aborted CME was reported on 
Feb. 12, but no flare was observed. The vector magnetic 
field was recorded b y the MEES/IVM dMickev et all 19961: 
LaBon te et al. I fl999l The high values of the current den- 
sity imply strongly sheared and twisted flu x bundles (see 
Regnier et al.l2002[|Regnier & Am ari 2004). Due to the ex- 
istence of highly twisted flux tubes (with more than 1 turn) 
and the stability of the reconstructed filament and sigmoid 
(with less than 1 turn), the authors concluded that the erup- 
tive phenomena was most likely to be due to the develop- 
ment of a kink instability in the highly twisted flux bundles; 
AR8210: observed on May 1, 1998 from 17:00 to 21:30 
UT, this is a newly emerged active region with a com - 
plex topology as described in iRegnier & Canfield (2006). 
An Ml. 2 flare was recorded on May 1, 1998 at 22:30 UT. 
The selected vector magnetogram (MEES/IVM) at 19:40 
UT was obser ved during a "quiet" period between two C- 
class flares. In Regnier & Canfield d2006h . the authors de- 
scribed the magnetic reconnection processes occurring dur- 
ing this time period and leading to a local reorganisation 
of the magnetic field. The reconnection processes are re- 
lated to the slow clockwise rotation of the main sunspot or 
a fast moving, newly emerged polarity. Following the time 
evolution during 4 hours, the authors showed that the free 
magnetic energy decreases during the flare over a period 
of about 15 min, and the total magnetic energy is slightly 
increased during this time period; 

AR9077: this corre sponds to the famous Bastille day 
flare in 2000 (e.g.. |Liu & Zhang||200U lYan et al.1 [2001 



Fletcher et alj 120011) . The vector magnetogram was 



recorded at 16:33 UT after the X5.7 flare which occurred at 
10:30 UT. The active region was still in the magnetic reor- 
ganisation phase after the flare and "post"-flare loops were 
observed in 195A TRACE EUV images. The flare was also 
associated with a CME; 

AR10486: this active region is responsible for the main 
eruptions observed during the Halloween events (26 Oct. 
to 4 Nov. 2003). The MEES/IVM vector magnetogram was 
recorded on October 27, 2003 at 18:36 UT before the X17.2 
flare which occurred at 11:10 UT on October 28. The flar- 
ing activity of this active region and the ass ociated CMEs 



have been extensively studied. For instance, iMetcalf et al 



(2005) have shown that the large magnetic energy budget 
(~ 3 10 33 erg) on Oct. 29 is enough to power the extreme 
activity of this active region. 



For these particular active regions, the reduction of the 
full Stokes vector to derive the magneti c field has already 
been detailed in several articles (e.g., IRegnier et al.1 120021 ; 
Regnier & Canfield! 120061) - the 180-degree ambiguity in the 
transverse component was solv ed by usin g the algorithm de - 
veloped in Canfiel d~etaD d 19931) (see also lMetcalf et al.ll2006h . 



2.2. Nonlinear force-free field 



- AR8151: observed on February 11, 1998 at 17:36 UT, Under coronal conditions, the magnetic force dominates the 
this is an old decaying active region (decreasing magnetic pressure gradient and gravity, and so we regard the corona 
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above active regions as being well described by th e force- 
free approximation (see e. g., recent reviews by iRegniei 
2007; IWiegelm ann 2008). Throughout this article, the coro- 
nal magnetic configurations are computed from the nonlin- 
ear force-free (nlff) approximation based on a vector poten- 
tial Grad-Rubin d 19581) m ethod by using the XTRAPOL code 
dAmariet al.ll 199711 1999b . The nlff field is governed by the fol- 
lowing equations: 



V A B = aB, 
B ■ Vff = 0, 



V ■ B = 0, 



(5) 
(6) 
(7) 



where B is the magnetic field vector in the domain Q. above 
the photosphere, £Q, and a is a function of space defined as the 
ratio of the vertical current density, J z and the vertical magnetic 
field component, B z : 



a 



1 

zT 



(dBy 

\ dx 



dB x 
8y 



(8) 



From Eqn. (O, a is constant along a field line. In terms of the 
magnetic field B, the Grad-Rubin iterative scheme can be writ- 
ten in two steps as follows. Firstly, 



B <fl) . Vq ,(») = o i n o 



(9) 
(10) 



where 60^ is defined as the domain on the photosphere for 
which B z is positive (+) or negative (-). This step corresponds 
to the transport of a along field lines, and secondly 



VAfi^'^aV in a 
V ■ B (n+l) =0 in O, 
Bf + %n = b z>0 , 



(11) 
(12) 
(13) 



ev en if the pho t ospher e is not ex actly force-free ac c ording 



to iMetcalf et al.1 d!995l) . Recently, IWiegelmann et al.l (|2006) 
have developed a pre-processing method in order to minimise 
the effects of non-force-f reeness and to obtain a more chro- 
mospheric magnetogram ( Wiegel mann et al. I l2008l) . This way 
of improving nlff field extrapolations is still under investiga- 
tion. Our nlff extrapolations have been successf ully compared 



to EUV, soft X-rays and visible observations ( Regn 



2002; 



Regnier & AmarJ f2004; 



Regnier & Canfield 



er et al. 

120061: 



Regnier & Fleck 2004; Regnier et al. 2005). 



For the sake of comparison, we also compute the potential 
field with the same boundary conditions. The potential field 
corresponds to the minimum energy state that a magnetic con- 
figuration can reach if all the currents are dissipated and while 
the distribution of the vertical field is held fixed. 

2.3. Coronal Alfven and sound speeds in an 
isothermal atmosphere 

In order to compute the coronal Alfven speeds, we first need 
to define the thermodynamic parameters such as temperature, 
density and pressure. As a first approximation, the corona is as- 
sumed to be an isothermal atmosphere satisfying a hydrostatic 
equilibrium: 



Vp+pg = 0. 



(16) 



In agreement with the Harvard-Smithonian model of the solar 
atmosphere, we suppose the assumption is reasonable above a 
height (zo) of about 5 Mm above the photosphere. Let us then 
compare the effect of a constant gravity (go) with a gravita- 
tional field (g(z)) decreasing with distance from the Sun. In an 
isothermal atmosphere, the pressure and density are either: 



p(z) = p exp(-^), and 

p(z) = po exp (-^) for g = go, 



(17) 



lim \B\ = 0, 



(14) 



which is the step required to update the magnetic configuration. 
The boundary conditions on the photosphere are given by the 
distribution b Zi o of B z on SQ. (see Eqn. (TT3l ) and by the distri- 
bution ao of a on 6Q. for a given polarity (see Eqn. ([Toll). We 
also impose that 



B n = 



(15) 



where S is the surface of the computational box, n refers to 
the normal component to the surface. These conditions imply 
that no field line can enter or leave the computational box, or 
in other words the active region being studied is magnetically 
isolated. 

We ensure that the total unsigned magnetic flux is bal- 
anced within the field-of-view by surrounding the vector mag- 
netograms by SOHO/MDI line-of-sight field. The smooth- 
ing function is applied in order to smooth out the transi- 
tion between the vector field and the MDI magnetogram. 
From the Grad-Rubin method, we obtain a nlff equilibrium 



P(z) = Po 
p(z) = po 

where H = 
1.38 10~ 23 J-KT 1 , ft 



ex p( H (jfe+zo) (« «))' 
ex p(~tf (rI+z ) (fe^)) 



and 



for g = g(z) 



So R l 



(18) 



(Ro+z> 



IcbT Kpnipgo) is the pressure scale-height (&g = 
0.6 for a fully ionized coronal plasma, 
nip = 1.67 10~ 27 kg and go = g(R&) = 274 m-s" 2 ), p and 
po are characteristic values of the pressure and density at zo- 
The pressure scale-height and the density po are the two free 
parameters of the model. Typical values are: H — 50 Mm for a 
temperature of 1 MK and po = 10 9 cirT 3 . 

In these models, the Alfven and sound speeds are given by: 



B(x,y,z) 

v A (x,y,z) = — ; exp 

V/"opo 

for a constant gravitational field and 



l z-zo \ 
\ 2H r 



B(x,y,z) 
v A (x,y,z) = — ; exp 



2H (R + zo) \R + z 



- ZO 



(19) 



(20) 



Regnier et al.: Coronal Alfven speeds in an isothermal atmosphere I. Global properties 



5 



(21) 



for a gravitational field varying with height, and 

c 2 = 2j^£ = M. 

PO pM p ' 

where B is the magnetic field strength and y — 1 in an isother- 
mal atmosphere. Therefore, the sound speed is constant in the 
corona and the Alfven speed can be defined locally knowing 
the magnetic field strength at each point of the coronal volume. 
Assuming that at infinity the magnetic field strength decays as 
1 jz 1 for a dipolar field, we notice that 

lim VA(x,y,z) - +oo, (22) 

for a constant gravitational field, while for a gravitational field 
decreasing with height, we obtain 



lim v A (x,y,z) = 0. 



(23) 



The assumption of a gravitational field varying with height is 
therefore more reasonable at high altitudes in order to avoid 
the Alfven speed becoming unbounded and to have a smooth 
transition with the solar wind. 

3. Study of a bipolar field 

3.1. Nonlinear force-free bipolar field 

Before applying this model to examples of solar active regions, 
we perform a parameter study using a nonlinear force-free 
bipolar field. The free parameters are the pressure scale-height 
(or temperature) and the density at the base of the corona. 





(b) 





■ 



Fig. 1. (a) Vertical magnetic field on the photosphere, (b) verti- 
cal current density distribution imposed in the positive polarity. 

The bipolar vertical magnetic field component is simply 
given by a Gaussian distribution with the same full-width at 
half-maximum for both polarities of 15 Mm (see Fig. QJ). The 
maximum field strength is set to 2000 G. The considered field- 
of-view has dimensions 150 Mmxl50 Mm assuming a square 
pixel of 1 Mm 2 . The separation between the two polarities is 
50 Mm. The magnetic flux is balanced and the total unsigned 
magnetic flux is 7.1 10 21 Mx. A typical bipolar magnetic field 
used in this study is depicted in Fig. [T^. The field-of-view is 
large enough to have no magnetic flux on the edges of the sur- 
face. We impose a ring distribution of current at the positive 
source of the form: 

J2- ' 



where r is measured from the centre of the source and Co is 
a constant which ensure a zero net current. We always take 
care to (i) keep a zero net current, and (ii) have enough pix- 
els to describe the steep gradient between the positive and neg- 
ative currents. A typical distribution used in this experiment 
is shown in Fig. QJ). We set J z q to 10 mA-m -2 . To satisfy the 
Grad-Rubin boundary-value problem, the current distribution 
is defined in only one polarity. This equilibrium corresponds to 
a single twisted flux tube with an excess magnetic energy of 
3% above the potential field. 



3.2. Magnetic field strength 

We first study the distribution of magnetic field strength with 
height for potential and nlff fields. We plot the average field 
strength as a function of height in Fig. [2] The average strength 
values start from 80 G at the photosphere and tend to zero. The 
rapid decrease of the field strength is due to the confinement 
of the twisted flux tube at a height lower than 30 Mm due to a 
separation between polarities of 50 Mm. There is a small de- 
parture of the nlff field values from the potential values due to 
the storage of magnetic ener gy in the corona associated with 
the twisted flux tube (see e.g, Regnier & Priest 2007b). 



100 



to 

IS 




J z (r) = 2 J z0 [r z - C ] exp 



(24) 



60 80 
Height (Mm) 

Fig. 2. Average field strength as a function of height (Mm) for 
the nlff field (solid line) and the potential field (dashed line). 
The grey area is enclosed between the minimum and maximum 
field strengths at a given height. The log-log plot is drawn in the 
top-right corner. The dot-dashed line indicates the z _1 decay of 
the field. 



In Fig. 12 we also plot the log-log distribution of the average 
field strength with height in order to determine the shape of the 
decaying magnetic field. For both the potential and nlff fields, 
the average magnetic field strength decays rapidly from the 
photosphere up to about 20 Mm and then the rate of decline 
of the field strength reduces with a slope index greater than -1 
as indicated by the dot-dashed line in the log-log plot of Fig. [2] 
We focus our interest on an index of -1 which corresponds to 
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the typical slope of a dipolar field average on a surface. The 
whole curve is well fitted by a function of the following form: 



B(z) = Bo exp 



ln 2 |- 

Uo 



-*(!P 



(25) 



as deduced from the log-log plot, where B is the average field 
streng th and z the height (Bo and zo being constants). iDere 
(11996 ) found that the magnetic field decay below 40 Mm is 



well described by an exponential decay. We will discuss these 
typical curves further in Section 

3.3. Alfven speed 



6000 




40 60 80 100 
Height (Mm) 



120 140 



Fig. 3. Average Alfven speed (km-s" 1 ) as a function of height 
(Mm) for H = 50 Mm and p = 10 9 cm 4 . Blue and red solid 
curves are from the nlff field for the constant and varying grav- 
ity models respectively. Blue and red dashed curves are from 
the potential field for the constant and varying gravity models 
respectively. The red and blue curves are very similar and al- 
most superimposed. The dark and light grey areas indicate the 
spread of Alfven speed values at a given height for the constant 
and varying gravity models respectively: only above 60 Mm, 
one can slightly notice the difference between the two gravity 
models. 

We derive the Alfven speeds from the model described 
above with H = 50 Mm (1 MK) andpo = 10 9 cirT 3 . In Fig. [3] 
we plot the average Alfven speed at a given height from both 
the potential field (dashed line) and the nonlinear force-free 
field (solid line). The Alfven speed values are valid only above 
5 Mm (vertical solid straight line). Under the force-free as- 
sumption, the plasma B is considered to be much less than 1 . 
The plasma B can be expressed in terms of the Alfven and 
sound speeds: Ic^lyv 2 ^ In an isothermal atmosphere, we then 
obtain a minimum value for the Alfven speed: 



2c 



(26) 



As the sound speed is 120 km-s" 1 for H — 50 Mm, we plot 
VA.min = 170 km-s" 1 in Fig.[3]and Fig. [5^ (dashed straight line). 
The average Alfven speed is similar for both gravity models in 
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Fig. 4. Dynamic range of the Alfven speed at 10 Mm (top row) 
and 60 Mm (bottom row) for two sets of free parameters (po, 
H) in (cm" 3 , Mm): (10 9 , 25), (10 8 , 50) from left to right. 

the case of the bipolar field. The curves are exactly the same up 
to a height z — H. Thus, we only report here the values for the 
constant gravity model. The average Alfven speed decreases 
with height suggesting that the magnetic field decreases faster 
than the square root of the density with height. The Alfven 
speed values range between 2400 km-s -1 at 10 Mm and 300 
km-s -1 at 140 Mm. The scatter of the Alfven speed (grey areas 
in Fig. [3]) strongly depends on the height: a maximum value 
of 24000 km-s" 1 at 10 Mm with a standard deviation of 3800 
km-s" 1 , and a maximum value of 990 km-s" 1 at 60 Mm with 
a standard deviation of 200 km-s" 1 (see Fig. By drawing 
the dynamic range of the Alfven speed at different heights (see 
Fig.|4]i, we show that the Alfven speed is larger where the field 
strength is larger whatever the density or the pressure scale- 
height: the two opposite polarities are clearly seen at 10 Mm 
and the loop top is seen at 60 Mm. 

We now study the effects of two free parameters: the den- 
sity at the base of the corona po (see Fig. |5^) and the pressure 
scale-height H (see Fig.[3J)). In Fig.|5h, we plot the evolution of 
average Alfven speed with height for several values of the den- 
sity po varying from 10 8 to 5 10 9 cm" 3 . As the magnetic field 
is the same, the average Alfven speed follows the density vari- 
ations which can be seen from the values of the Alfven speed 
at 5 Mm: 1560 km-s" 1 for p = 5 10 9 cm" 3 , 3500 km-s" 1 for 
po = 10 9 cm" 3 , 4950 km-s" 1 forp = 5 10 8 cm" 3 and 11070 
km-s" 1 for po = 10 8 cm -3 . We note that there is a minimum of 
th ese curves a t about 135 Mm. This minimum was also noticed 
bvlLin1(l2002) but at about 0.5 R (~ 350 Mm) for an isothermal 



atmosphere. From Eqn. ( TT9b . the minimum corresponds to the 
location where the variations of the gas pressure are becoming 
more important than those of the magnetic field (that does not 
mean that the plasma B becomes greater than 1 !). We conclude 
that this minimum is strongly related to the distribution of the 
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Fig. 5. Average Alfven speed (km-s -1 ) as a function of height 
(Mm) for the nonlinear force-free field (solid line) and the 
potential field (dashed line), (a): for H — 50 Mm and po = 
[5 10 9 , 10 9 ,5 10 8 , 10 8 ] cm" 3 (from left to right); (b): for H = 
[100, 50, 25] Mm (from left to right) and p = 10 9 cm 3 . Red 
(resp. blue) curves are for the varying (resp. constant) gravity 
model. 

magnetic field: from Eqn. [19] we deduce that the minimum of 
the curves is obtained for a constant gravitational field when 
the magnetic field satisfies 

dB _ B 

~dz ~ ~2h' 



(27) 



which depends only on the pressure scale-height. For a gravi- 
tational field varying with height, the minimum occurs where 



B 



dB 

~dz 2H\R Q +z 



(28) 



We can already mention that Eqn. J27b is similar to the 
equation derived by iFilippov & Denl (I2000f) for the height of 
unstable prominences. In Fig. |5j), we plot the average Alfven 
speed with height for several values of the pressure scale-height 
varying from 25 to 100 Mm (from 0.5 to 2 MK, resp.) giving a 
sound speed between 83 km-s -1 and 165 km-s -1 , respectively. 
Therefore, the Alfven speed curve for H = 100 Mm is less than 
VA,min above 100 Mm, and the curve for H = 25 Mm is way 



above whatever the height. The minimum of the curves 
is decreasing with height when H decreases. The Alfven speed 
values increase significantly when H decreases: at 100 Mm, the 
Alfven speed is 200 km s 1 for H = 100 Mm, 400 km s 1 for 
H = 50 Mm and 1000 km s 1 for H = 25 Mm. 

From this study, we conclude that a density at the base of 
the corona of 10 9 cm -3 and a pressure scale-height of 50 Mm 
are a reasonable choice of parameters to model an isothermal 
corona. We will adopt these values in the following Sections. 



4. Active regions 

4.1. Magnetic field strength in active regions 

Before computing the Alfven speed, we study the changes of 
magnetic field strength depending on the nature of the observed 
active region and on the extrapolation model. We first compute 
the potential and nonlinear force-free magnetic fields associ- 
ated with the four active regions described in Section l2~T1 Then 
we derive the magnetic field strength average in an xy-plane as 
a function of height. In Fig. [6] we plot the average field strength 
for the potential field (dashed line) and the nonlinear force-free 
field (solid line). The grey area is the variation of field strength 
for the nonlinear force-free field at a given height. 

For AR 8151 (Fig.|6]top left), the average nlff field strength 
ranges from 55 G on the photosphere to 9 G at 170 Mm. 
Assuming that the decrease of the field strength should be ex- 
ponential, there is an excess in the values of the field strength 
between 40 Mm an d 100 Mm where the t wisted flux bun- 
dles ar e located (see Regnie r & Amaril[2004 ; Regnier & Priestl 
2007bh . 

For AR 8210 (Fig. top right), the average nlff field 
strength varies from 45 G on the photosphere to small values 
above 300 Mm. The nlff curve (solid line) follows an exponen- 
tial decay as for the potential field (dashed line). It is noticeable 
that the scatter of the field strength below 50 Mm is of about 
two orders of magnitude. Even i f the storage of magnetic en- 
ergy is focussed near the surface (Regnier & Priest 2007b), the 
largest departure from potential occurs between 50 Mm and 
300 Mm. Below 50 Mm, the average values of the nlff field 
strength are similar to the potential ones but the spatial distri- 
bution is different: the free magnetic energy above the potential 
field is only 2.5% of the total magnetic energy and the energy is 
st ored near topolog i cal elem ents in the nlff field, as mentioned 
in 



Regnier & Priest d 007alTbh . 
For AR 9077 (Fig. |6] bottom left), the total magnetic flux 
on the photosphere is more important than for the above ex- 
amples. The average field strength varies from 300 G on the 
photosphere to 90 G at 110 Mm. The average strength from 
the nlff field follows the same decay as for the potential field. 
The active region was observed after a X5.7 flare and exhibits 
post-flare loops. The variation of the field strength with height 
indicates very little storage of magnetic energy in the corona 
due to twisted flux tubes contrary to AR 8151. 

For AR 10486 (Fig. |6] bottom right), the average field 
strength varies from 250 G on the photosphere to 5 G at 250 
Mm. It is noticeable that the magnetic energy can be stored in 
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Fig. 6. Average magnetic field strength as a function of height for the potential field (dashed line) and for the nonlinear force-free 
field (solid line). The grey area ranges from the minimum to the maximum field strength of the nonlinear force-free field. The 
log-log plots are drawn in the top-right corners. The dot-dashed line indicates the z _1 decay of the field. 



the corona up to 150 Mm as seen by the departure from the 
potential field (dashed curve). 

By studying the global properties of the magnetic field 
strength above several active regions, we conclude that: 

(i) the main magnetic flux concentration occurs near the pho- 
tosphere up to a height of 50 Mm; 

(ii) a significant amount of magnetic energy can be stored in the 
corona (between 50 and 150 Mm) in an active region with 
highly twisted flux bundles (AR 8151) or with a complex 
topology associated with a strong photospheric flux (AR 
10486); 

(iii) the potential and nlff fields have different behaviours de- 
pending on the nature of the field: strong departure from 
the potential field for active regions with highly twisted flux 
tubes or complex topology associated with a high activity 
level, and small departure for active regions after a strong 
flare or with a complex topology not related to important 
flare activity. 



As pointed out by iKoutchmy (2004), there are few mea- 
surements of the magnetic field strength in the solar atmo- 
sphere. The above study shows that the average values at a 



given height strongly depend on the nature of the active re- 
gion (the total magnetic flux and the distribution of polarities 
on the photosphere) and on the nature of the coronal magnetic 
field (potential or nlff fields). In Section [5] we will discuss the 
implications for different solar models. 

Regarding the relationship between eruption and magnetic 
field decay, the log-log plots (see Fig. [6]) show several interest- 
ing features: 



(i) the z _I decay line is tangent to the curve of the nlff field 
at about 80 Mm for both active regions (AR 8151 and AR 
10486) associated with an eruption causing a large-scale 
magnetic field reorganisation, whilst for a post-eruptive ac- 
tive region (AR 9077) and the one associated with confined 
flares (AR 8210), the z" 1 line is at about 20 Mm; 

(ii) the weak slope of the log-log plots below the z _1 line in- 
dicates the strong influenc e of the current systems as in 
van Tend & Kuperusl dl978l) . 



The implications for eruption models are developed in 
Section|5] 
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4.2. Alfven speed in solar active regions 

We apply the above models to the active regions described in 
Section [2T| To avoid confusion, we only report here the Alfven 
speed values for the constant gravity model. In Fig. [7] we plot 
the average Alfven speed values in the xy-plane as a function 
of height derived from the potential field (dashed curves) and 
nlff field (solid curves) for the different active regions. First, the 
average Alfven speed at the base of the corona varies from 2- 
3000 kms 1 for AR 8151 andAR8210to 13000 kms 1 forAR 
9077 and AR 10486. It is strongly related to the total magnetic 
flux of the active region on the photosphere. Above the coronal 
base, the variation of the Alfven speed with height is greatly 
influenced by the nature of the field: the topology, the geometry 
of field lines and the effects of currents. 

For AR 8151, the Alfven speeds in the nlff field can reach 
up to 52000 km-s -1 near the base of the corona (for an average 
value of 3000 km-s~'). Then the maximum Alfven speed val- 
ues decrease to 6000 km-s -1 at 50 Mm to finally reach the value 
of 5000 km-s -1 above 150 Mm. The evolution of the average 
values is almost constant at 3000 km-s~'. Nevertheless there is 
a minimum of Alfven s peed of 2500 km-s" 1 at 2 5 Mm. This 
re sult is consistent with IWarmuth & Mannl ( 12005b . According 



to lRegnier & Priest (2007b), the minimum is below the height 



where the twisted flux tubes are located and the magnetic en- 
ergy is deposited in this active region. The evolution of the av- 
erage speeds in the potential field (dashed curve) is strictly de- 
creasing from 2500 km-s -1 to 1500 km-s~'. Note that accord- 
ing to Section[3] the average speeds should increase above 170 
Mm when the density dominates. Both potential and nlff av- 
erage speed curves are above the /3 < 1 limit (dashed line). 
Even if the average values of the Alfven speed for the potential 
and nlff fields are not far from each other in the low corona, 
there is a factor of about 2 difference above 80 Mm. In Fig. [8] 
we draw the dynamic range of Alfven speed at two different 
heights (top: 10 Mm, bottom: 60 Mm) for both potential (left) 
and nlff (right) fields. At both heights and for both models, the 
greater is the magnetic field strength, the larger is the Alfven 
speed. Even if the dynamic range is similar for the potential 
and nlff fields at 10 Mm, the Alfven speed varies by a factor of 
2 at 60 Mm due to the presence of twisted flux bundles along 
the polarity inversion line. The average Alfven speed decreases 
more rapidly for the varying gravity model: at a height of 150 
Mm, the average Alfven speed is 2400 km-s -1 for a varying 
gravity and 3000 km-s -1 for a constant gravity. 

For AR 8210, the maximum Alfven speed at the base of 
the corona is about 70000 km s 1 . Even if the total mag- 
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Fig. 8. Alfven speed dynamic range in AR 8151 for the poten- 
tial field (left column) and the nonlinear force-free field (right 
column) at two different heights, 10 Mm (top) and 60 Mm (bot- 
tom). 
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average speed for the nlff field (potential resp.) is about 650 
km-s -1 (400km-s _I resp.) which is indeed of factor of 1.5 more 
for the nlff field. The dynamic range of the Alfven speed is de- 
picted in Fig. [9] There are very few changes between the dy- 
namic range images at a given height. We then suggest that the 
complex topology of AR 8210 does not significantly influence 
the Alfven speed distribution. This conclusion is also supported 
by the fact that the magnetic energy is deposited near the pho- 
tospheric boundary and not in the corona. Since the average 
Alfven speed decreases faster for the varying gravity model, 
the Alfven speed becomes less than the coronal sound speed at 
about 300 Mm, whilst the Alfven speed stays above the plasma 
/3 = 1 limit for the constant gravity model. 




Fig. 10. Alfven speed dynamic range in AR 9077 for the poten- 
tial field (left column) and the nonlinear force-free field (right 
column) at two different heights, 10 Mm (top) and 60 Mm (bot- 
tom). 



Fig. 9. Alfven speed dynamic range in AR 8210 for the poten- 
tial field (left column) and the nonlinear force-free field (right 
column) at two different heights, 10 Mm (top) and 60 Mm (bot- 
tom). 



netic flux is comparable to the flux of AR 8151, the maxi- 
mum field strength is larg er in the main sunspot of AR 8210 
(Reg nier & Canfieldl 120061) which explains the larger value of 
the maximum Alfven speed. The average speed for the nlff field 
decreases rapidly from the base of the corona to about 70 Mm 
from 2200 km- s~ 1 to 600 km- s~ 1 . The fast decrease is explained 
by the nature of the magnetic configuration: the confined mag- 
netic field in the active regio n, and the storage of mag netic 
energy in the low corona (see lRegnier & Canfieldl 120061) . The 
average speed for the potential field follows the same fast de- 
crease below 70 Mm before reaching a minimum at about 150 
Mm. Above this height, the average speed slowly increases 
when the density dominates the field strength. At 100 Mm, the 



For AR 9077, the total magnetic flux included in the field- 
of-view is larger than for the previous two examples. Therefore, 
even if the maximum Alfven speed at the base of the corona 
is moderate (~ 46000 km-s -1 ), the average speed is relatively 
large (13000 km-s -1 ). Both average curves for the potential 
and nlff field are rapidly and strictly decreasing towards 1500 
km-s -1 . Contrary to the previous example, there is no evidence 
of storage of energy in the middle corona. The strong flare and 
the CME associated with this active region are related to a com- 
plex topology and the mass supply is prov ided by a low fila- 
ment on the side of the active region (e.g., Zhang 2002). The 
reconstructed active region is performed from a magnetogram 
recorded after the flare. As suggested by Fig. [10] the dynamic 
ranges are quite different from one model to another, whatever 
the altitude we look at, even if the Alfven speed values are com- 
parable. Therefore the post-eruptive configuration is strongly 
influenced by the electric currents flowing along field lines de- 
spite a magnetic geometry which looks similar to potential. For 
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What can we learn from the estimate of the coronal Alfven 
speed for some key solar physics problems (coronal heat- 
ing, magnetic reconnection, wave propagation, initiation 
and development of CMEs)? 




Fig. 11. Alfven speed dynamic range in AR 10486 for the po- 
tential field (left column) and the nonlinear force-free field 
(right column) at two different heights, 10 Mm (top) and 60 
Mm (bottom). 



both gravity models, the average Alfven speed curves are sim- 
ilar indicating that, during this relaxation phase, gravity does 
not greatly influence the magnetic field and the plasma. 

For AR 10486, the maximum Alfven speed at the base of 
the corona is estimated to be 130000 km-s~', which is about 
one-third of the speed of light. The average speed for both po- 
tential and n/jf fields is about 13000 km-s~'. The average speed 
for the nlff field is decreasing rapidly from 13000 km-s -1 at 
5 Mm to 6000 km-s~' at 50 Mm and then reaches an almost 
constant value of 5000 km-s" 1 . In the view of the previous 
examples, we interpret this fact as a consequence of the stor- 
age of magnetic energy in the corona at a height of about 
50 Mm ei ther by twisted flux tub es or by sheared arcades. 
Following iRegnier & Fleckl (2004), the existence of sheared 
arcades in a complex topology is confirmed. The effects of 
the transverse field components and so of the electric currents 
are shown in Fig. QT| The location of strong magnetic field 
strength is different from one model to another showing the 
com plex distribution of the magnetic field on the photosphere 
(e.g. jRegnier & Fleckl2004l;lMandrini et alj2006l) . The Alfven 
speed for the constant gravity model is almost twice the value 
for the varying gravitational field high in the corona. 



5. Discussions and Conclusions 

In this paper, we address two questions related to the global 
properties of the magnetic field and the plasma in the corona: 

1- How does the magnetic field strength and Alfven speed 
vary with height in the solar corona? 



In order to tackle the above questions, we have developed 
a method to determine the plasma properties in a magnetic 
field above active regions. We assume that an active region can 
be modelled in two steps: a nonlinear force-free equilibrium 
for the magnetic field, and a hydrostatic equilibrium between 
pressure and gravity forces describing the plasma properties. 
To solve the latter equation, the corona is assumed here to be 
isothermal with either a constant gravity or a gravitational field 
varying with height. 



In Sections 13.21 and 14.11 we investigated the evolution of 
the average field strength with height. Interestingly, the log-log 
plots showed that the decay of the magnetic field with height 
does not follow a power law, even for a potential field. This 
fact suggests that the decay of the field strength is strongly 
influenced by the geometry of the field (e.g., toroidal geome- 
try for the bipolar field), by the distribution of electric currents 
and their typical height of influence, and by the complexity of 
the magnetic distribution on the photosphere for active regions. 
This is a well-know n property of the magnetic field. Based on 
the da ta reported by lPoletto et al.l (1 19751) . Ivan Tend & Kuperus 
(1978) plotted the measured magnetic field strength as a func- 
tion of height for three different active regions and assum- 
ing current- free magnetic configurations. They obtained simi- 
lar curves to those shown in this paper in which the log-log plot 
can be fitted by a parabola. In addition, they fitted the distribu- 
tion of field strength by two power laws of index -1 and -3 sug- 
gested by their theoretical work: - 1 for a line current along the 
inversion li ne of the magnetic field and -3 for a dipolar potential 
field decay, van Ten d & Kuperus dl978l) concluded that (i) the 
electric currents strongly influence the magnetic configuration 
when the power law index is less than -1, (ii) above a certain 
height the magnetic field decays rapidly (probably with a power 
law index of -3). More importantly, the authors concluded that, 
if the magnetic field below the z decay is decaying with a 
power law index less than -1, then the configuration is stable. 
This condition can be compare d with the conditions required to 
trigger a ki nk instability (e.g.. | Hood & Priesdl 1 98 lh or a torus 
instability dKliem & Torokll2006l) . Recently. iLiul d2008l) inves- 
tigated the evolution of the magnetic field strength with height 
and compared with the nature of the eruptive phenomenon. He 
derived the magnetic field from a large-scale potential field ex- 
trapolation ba sed on the PFSS m odel (Potential Field Source 
Surface, e.g JSchatten et al. 1969). For heights between 42 Mm 
and 140 Mm, the power law index is about -1.5 whatever the 
onset mechanism. Combined with our study, this result sug- 
gests that the behaviour of the magnetic field in the low at- 
mosphere plays the main role in triggering flares and CMEs, 
and also that electric currents flowing along field lines as in the 

model affect strongly the power law index and therefore the 
stability of the magnetic configuration. 



12 



Regnier et al.: Coronal Alfven speeds in an isothermal atmosphere I. Global properties 



100 r 



■a 



Co 



1 



Height (Mm) 



M,=0.1 



M =0.01 - 



M =0.001 



100.0 



■a 10.0 



Co 



1 .0 



0. 





AR8210 















100 200 300 

Height (Mm) 



■a 



1000 



100 r 



CO 




1000 



40 60 80 

Height (Mm) 



Co 
o 

s? 



100 




100 150 200 
Height (Mm) 



Fig. 12. Average inflow speeds as a function of height for the nonlinear force-free field (solid line) and potential field (dashed 
line), with a reconnection rate M e varying from 10~ 3 to 10 _1 . Red (resp. blue) curves are for the varying (resp. constant) gravity 
model. 



5.1. Reconnection rate 

As mentioned in SectionQ] the Alfven speed is an important pa- 
rameter in understanding the magnetic reconnection processes. 
The reconnection is characterised by the reconnection rate (or 
external Alfven Mach number) as follows: 



M e = —, 
va 



(29) 



where v, is the inflow speed. This definition does not depend 
on the nature of the diffusion region where the magnetic re- 
connection occurs. To distinguish between different regimes of 
reconnection, the magnetic Reynolds number is introduced as 
follows: 

V 



L e VA 



(30) 



where rj is the magnetic diffusivity, and L„ is a t ypical length of 
the diffusion region (e.g. JPriest & Forbesl2000l) . If Rz} c < M e < 



— 1/1 

R„J , we have a slow reconnection regime of the Sweet-Parker 
type. If M e > Rm l J 2 , the reconnection regime is fast such as in 



the Petschekl l 1964) model. For a typical astrophysical plasma, 
the reconnection rate in the fast regime is in the range 0.01 to 
0. 1 . The square of the reconnection rate is also the ratio of ki- 
netic energy to magnetic energy of the inflow region assuming 



that the diffusion region is small enough to have the same field 
strength in the inflow and outflow regions. For a given inflow 
speed, the higher is the magnetic energy stored in the region, 
the smaller is the reconnection rate. 

If we assume t hat the reconnection rate is between 10~ 3 
and 10~' (see e.g., iNarukage & Shibata 2006), we can infer 
the inflow speed required to trigger the reconnection. The re- 
sults for the four studied active regions are summarized in 
Fig. [12] where we plot the inflow speed as a function of 
height for different reconnection rates (M e - 0.1,0.01,0.001) 
for both potential (dashed curves) and nlff (solid curves) 
fields. The variation with height of the inflow speed is con- 
sistent with the variation of the Alfven speed as described 
in Fig. 13 From these curves, we can deduce a threshold 
on the reconnection rate in order to obtain inflow speeds in 
agreement with observations (e .g., Narukage & Shi batall2006l : 
Nagashi ma & Yokovamal200"6l) . For instance, if a reconnection 



process occurs in AR 8151, the reconnection rate should be of 
the order or below 0.01 which gives typical inflow speeds be- 
low 50 km-s -1 . 

We also note that the minimum of the inflow speed or 
Alfven speed (see Section [3]) is the location in the low corona 
where the reconnection process is more likely to happen. The 
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Fig. 13. Plasma fi for the nonlinear force-free field (solid line) and the potential field (dashed line) as a function of height for four 
active regions. Red (resp. blue) curves are for the varying (resp. constant) gravity model. 



conclusion is valid in a statistical sense because we only con- 
sider average values of the inflow speed and because, in or- 
der for it to occur, the reconnection process needs a particular 
magnetic topology and a diffusion region which are beyond the 
scope of this article. 

5.2. Coronal plasma f3 

An important quantity in plasma physics is the plasma fi, the ra- 
tio of the gas pressure to the magnetic pressure. When fi <sz 1, 
the plasma is dominated by the magnetic field. Nevertheless, 
the plasm a fi in the solar atmosphere varies greatly with height 
dGarvl200ll) . In addition, we note that the plasma fi varies from 
one structure to another and even coronal structures, such as fil- 
aments, can have fi > 1. As mentioned in Section[3] the plasma 
fi can be expressed as follows: 



2c2 
P=—t 



(31) 



for an isothermal atmosphere (y = 1). In Fig. [13] we plot the 
variation of the average plasma fi with height for the four ac- 
tive regions described in Section 12.11 For both potential and 
n/j(f models, the plasma fi is less than 1 . This result is consistent 



with the typical nature of the coronal active region magnetic 
field. The plasma fi is larger for the potential field because the 
average Alfven speed decreases with height more rapidly than 
for the nlff field. The maximum of the plasma fi is located at the 
same height as the minimum of Alfven speed (see Section [3~3l 
and also satisfies Eqn. (1271 1. The plasma fi tends to zero for 
a constant gravitational field, whilst it increases towards 1 for 
the varying gravity model. The latte r is then in ag reement with 
plasma fi measurements reported by iGarvl (12001). 

5.3. Average Alfven transit time 

The Alfven transit time is an important quantity for the prop- 
agation of Alfven waves and for the relaxation of a mag- 
netic configuration. In this study of global properties of Alfven 
speeds, we derive the average Alfven transit time associated 
with two different lengths: 



and 



'a 



// 

VA 



nL 



V/l 



(32) 



(33) 
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Fig. 14. Average Alfven transit time for the potential (dashed line) and nonlinear force-free (solid line) fields using as typical scale 
the pressure scale-height H (bottom curves) or the length of a semi-circular loop with a diameter of the scale of the computational 
box (top curves). Red (resp. blue) curves are for the varying (resp. constant) gravity model. 



where v& is the Alfven speed, H is the pressure scale-height 
and L is the length of the computational box. We suppose 
(resp. t£ ) is a minimum (resp. maximum) of the Alfven tran- 
sit time. The larger the Alfven transit time is, the more stable 
the equilibrium is. In Fig. [14] we plot the different Alfven tran- 
sit times for the four active regions. At the base of the corona 
(5 Mm), the maximum Alfven time is 750s for AR 8 1 5 1 , 3600s 
for AR 8210, 70s for AR 9077, 300s for AR 10486. These re- 
sults mean that the most stable active region is AR 8210 associ- 
ated with confined C-class flares, and the least stable one is AR 
9077 exhibiting post-flare loops. It is worth noticing that the av- 
erage Alfven transit time is an interesting quantity only when 
discussing the global or statistical properties of magnetic con- 
figurations. A more appropriate quantity is the Alfven transit 
time along a single field line which is beyond the scope of this 
article. Of importance for MHD modelling, the Alfven transit 
time along field lines is not a constant with height which makes 
it difficult to estimate the scaling of the time evolution with this 
quantity. 

Despite the simple assumptions of this model, we have de- 
rived interesting properties of the Alfven speed in the solar 



corona. In a forthcoming paper, we will focus on the proper- 
ties of individual field lines. 
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